% clc;
% clear;

load('CE_EDC.mat');
CE = EDC;
load('AE_EDC_2.mat');
AE = EDC;
eddyArray = EddyDatasetClass.eddyArrayAddByDict([],CE.Dict,1);
eddyArray = EddyDatasetClass.eddyArrayAddByDict(eddyArray,AE.Dict,-1);
t = datetime(2011,8,5);
eddyArray = eddyArray(find([eddyArray.time] == t));
% load('ETTEC.mat');
extent = [0 25;108 122]; 
EDrawC = EddyDrawClass('extent',extent);
%18 50m
%23 100m
%25 150m
%27 200m
EDrawC.drawEddyWithSSTA(eddyArray, 33);

surfece = ssta(:,:,2);
%%
% % 提取涡旋数组
% eddyArray = EddyDatasetClass.eddyArrayAddByTraj([], EDC.traj, -1);
% 
% % 取测试数组
% eddys = eddyArray(10000:11000);
% 
% EAMSC = EddyArrayMatchSSTAClass(eddys);
% 
% % 求涡旋平均三维结构和平均半径
% [SSTA, R] = EAMSC.mean_eddyMatchedSSTA_Radis();
% 
% % 画涡旋剖面
% EDrawC.DepthProfileDraw(SSTA);
% % 画涡旋截面
% EDrawC.CrossSectionDraw(SSTA);

%%
% 设定区域必要画图条件
% EDrawC.extentSet(ETTEC.extent);
% EDrawC.intervalSet(ETTEC.interval);

% EDrawC.extentSet(ETTEC.extent);
% EDrawC.intervalSet(ETTEC.interval);
% 
% area_mean = ETTEC.areaThermalmeasure();
% [MHT, ZHT] = ETTEC.eddyHeatTransportMeasure();
% EDrawC.ThermalSectionDraw(area_mean,'depth',10,'reshapeInterval',0.25);
% EDrawC.EHTDraw(MHT, ZHT);
% save('AE_EddyDraw_2_1X1.mat');